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We study the magnetic soliton dynamics of spinor Bose-Einstein condensates in 
an optical lattice which results in an effective Hamiltonian of anisotropic pseudospin 
chain. An equation of nonlinear Schrodinger type is derived and exact magnetic 
soliton solutions are obtained analytically by means of Hirota method. Our results 
show that the critical external field is needed for creating the magnetic soliton in 
spinor Bose-Einstein condensates. The soliton size, velocity and shape frequency can 
be controlled in practical experiment by adjusting the magnetic field. Moreover, the 
elastic collision of two solitons is investigated in detail. 
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I. INTRODUCTION 



The realization of spinor Bose-Einstein condensates (BECs) opens a useful tool to under- 
stand and to confirm the dynamics of periodic structures in solid state physics. Recently, 
the research of spinor BECs trapped in optical potentials have received much attention both 
experimentally jl], O] and theoretically jsj. Due to the internal degrees of freedom for the 
hyperfine spin of the atoms, spinor BECs bring forth a rich variety of phenomena such as 
spin domains 0, and textures [si]. When the potential valley is so deep that the individual 
sites are mutually independent, spinor BECs at each lattice site behave like spin magnets 
and can interact with each other through both the light-induced and the static, magnetic 
dipole-dipole interactions. These site-to-site dipolar interactions can cause the ferromagnetic 
phase transition 7, sj leading to a "macroscopic" ma gne tization of the condensate array, the 



spin- wave like excitation |7HlO| and magnetic soliton [Hi, 112 1 analogous to the spin- wave and 
magnetic soliton in a ferromagnetic spin chain. For a practical spin chain, the site-to-site 
interaction is caused mainly by the exchange interaction, while the dipole-dipole interaction 
is negligibly small. For the spinor BECs in the optical lattice, the exchange interaction is 



2 



absent. The individual spin magnets are coupled by the magnetic and the light-induced 
dipole-dipole interactions jo], Q which are no longer negligible due to the large number of 
atoms N at each lattice site, typically of the order of 1000 or more. Therefore, the spinor 
BECs in an optical lattice offer a totally new environment to study spin dynamics in periodic 
structures. The magnetic soliton excited by the interaction between the spin waves is an 
important and interesting phenomenon in spinor BECs. The Heisenberg model of spin-spin 
interactions plays a significant role in understanding many complex magnetic structures in 
solids. It explains the existence of ferromagnetism and antiferromagnetism at temperatures 
below the Curie temperature. The magnetic soliton which describes localized magne- 



tization, is an important nonlinear excitation in the Heisenberg model |14j-ll7l|. However, 
the generation of controllable solitons is an extremely difficult task due to the complexity of 
the conventional magnetic materials. The spinor BECs seems an ideal system to serve as a 
new test ground for studying the nonlinear excitations of spin waves both theoretically and 
experimentally. 

In this paper, we demonstrate that the magnetic soliton and elastic soliton collision are 
admitted for spinor BECs in a one-dimensional optical lattice. By means of Hirota method 
we obtain the analytical dark soliton solutions, and also discuss the soliton interaction in 
detail. The outline of this paper is organized as follows: In Sec. II the Solitons in a Spin 
Chain of Atomic BECs is investigated in detail. Next, we obtain the one-soliton solution of 
spinor BECs in an optical lattice. In Sec. VI, the general two-soliton solution is obtained. 
Analysis reveals that elastic soliton collision occurs and there is a phase exchange during 
collision. Finally, our concluding remarks are given in Sec. V. 



II. SOLITONS IN A SPIN CHAIN OF ATOMIC BECS 



The Hamiltonian describing an F = 1 spinor condensate trapped in an optical lattice, 
which is subject to the magnetic dipole-dipole interactions and is coupled to an external 
magnetic field via the magnetic dipole Hamiltonian Hb has the form as |3|, |j, [g, [3] 



H = Hq + Hd + Hb, 



(1) 
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where the first term Hq describes the interaction of the atoms with the lattice potential 
Ul (r) and the ground-state collisions which can be written as 

Ho = J2f ^rV^lW[-^ + f/L(r)]^Jr) 

where tp^ (r) is the field annihilation operator for an atom in the hyperfine state 
1/ = l,mf = a), the indices z/ = 0, ±1 denote the Zeeman sublevels of the ground 

state, the parameters and characterize the short-range spin-independent and spin- 
changing s-wave collisions, respectively. When the optical lattice potential is deep enough 
there is no spatial overlap between the condensates at different lattice sites. Under this 
condition the tight-binding approximation is reasonable that the atomic field operator can 
be expanded as -ip (r) = J2n X]a=o±i (n) 0„ (r), where n labels the lattice sites, 0„(r) is 
the Hartree wave function of the condensate for the nth microtrap and the operators aa{n) 
satisfy the bosonic commutation relations [aQ,(n), dj3(/)] = Sa/sSni- It is assumed that all 
Zeeman components share the same spatial wave function. If the condensates at each lattice 
site contain the same number of atoms A^, the ground-state wave functions for different sites 
have the same form 0„ (r) = / (r — r„). In this case the dipole-dipole interaction potential 
Hd in Eq. ([1]) can be expressed by 

Hd = ^Yl |3 ll^i ■ T^i - 3 Cti ■ fy) itj ■ fii)] , 

where the parameter jti = 78^ is the magnetic dipole moment at site i, with Sj = 
al^{i)Fai3ai3{i) being the total angular momentum operator and 7 the gyromagnetic ratio, 
Ti denotes the coordinate of the ith site, r^j = — r^, fjj = Vij/ |rjj|, /Xq is the vacuum 
permeability. Also the Zeeman energy Hb can be described by 

i 

where B is the external magnetic field. 

In this paper we consider a one-dimensional optical lattice along the z-direction, which 
is also chosen as the quantization axis. In the absence of spatial overlap between individual 



4 



condensates, and neglecting unimportant constants, the effective spin Hamiltonian can be 
71, y as 



constructed 



H 



(2) 



which determine the dynamics of the magnetic dipole moment of spinor BECs in an optical 
lattice, where the parameter = (l/2)Aa/ rf'^r |0„(r)|'', here \a characterizes the spin- 
dependent s-wave collisions, and Jij = 7^/io/(47r |rjj|^). The first term in Eq. (|2]) represent 
the spin-dependent interatomic collisions at a given site, the second and the third terms 
describe the site-to-site spin coupling induced by the static magnetic field dipole-dipole 
interaction. Without this interaction terms, the ground state of the Hamiltonian is 
1^5) = |A^, — A^) where = Ni is the total atomic number in the lattice. The total spin 
at site i has the expectation value (^S^^ = —Nih. Due to the large factor Ni, the magnetic 
dipole-dipole interaction in the optical lattice can not be neglected. After the site-to-site 
coupling is considered, the transfer of the transverse spin excitation from site to site is 
allowed, resulting in the distortion of the ground state spin structure. This distortion can 
propagate and hence generate magnetic soliton or spin wave along the atomic spin chain. 

reported earlier is 



It should be noted that the Holstein-Primakoff transformations [18 1 
an useful tool to clear the physical meaning of spin waves and magnetic soliton. For this 
transformation a local spin-deviation operator is introduced firstly, i.e., fi = S — with 
the eigenvalues n = S — m in which it can be seen that the increasing m decreases n 
and vice versa. For the state |n), the \n) = a/ (25* — n){l + n) \n + 1) and \n) = 



{2S — n) {n — 1) \n — 1), it can define the creation and annihilation operators a and 
which satisfy the following boson commutator relation [a, a+] = 1, [a, a] = [a+, a+] = 0. The 



result of a and a"*" operate on the state \n) is a \n) = ^Jn\n — 1), a+ |n) = -y/n + 1 |n + 1), 
and a^a \n) = n\n) . With the above relations, the spin operators S can be expressed by the 
Bose operators a"*" and a, i.e., Holstein-Primakoff transformations 18| as follows 



5+ 



[\J2S — a+a)a, 
a+{V2S-a+d) 
{S — a^d). 



(3) 



In fact the spin waves and magnetic soliton show the small distortion of the ground-state 



5 



spin structure, i.e., n S. In this case we can expand the Eq. as 

S- = a+v^(l-^ + ...), 

S, = (5-a+a). (4) 

Substituting the Eq. (jl]) into the Hamiltonian ([2]) and keeping terms through fourth order 
in d and d^, we can get 



djdj + d|dj 



H = NX'^S (S + l)- -fNSB, + 75, ^ d|di + ^ J] 5Z 

— - Jij ^djdjdjdj + djdjdjdj + d|djdjdj + d|djdidj^ 

—2 Jij i^S^ — S ^djdj + d|dj j + d|djd]dj j . (5) 

Under the spin coherent state and using the time-dependent variation principle, the nonhnear 
operator motion equation of Hamiltonian ([5]) can be transformed into an equation for the 
probability amplitude V'fe = {'^\0'k\'^) which describes the nonlinear dynamics of coherent 
spin excitations on the lattice k. 
d_ 



Yl i'^^li^ki^k+s + i^l+si^k+si^k+s + i^ki^ki^l+s) (6) 

^ s 

where 5 = ±1, i.e., we consider only the nearest-neighbor interactions which is a good 
approximation [19] for the BECs in a one dimensional optical lattice as the large lattice 
constant. 



III. ONE SOLITON SOLUTION 

When the optical lattice is infinitely long and the spin excitations are in the long- 
wavelength limit, tp;^, ^ (^)^) the continuum limit approximation, we have 

- + + 8S^^ + 2S^,, - 12^ IV^I^ = 0. (7) 

In fact, Eq. ([7]) has dark soliton solutions which can be obtained by many methods. Here 
we take the Hirota bilinear transformation to get the exact one- and two-soliton solutions 
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of Eq. ([7]). For this method it apply the direct transformation to the nonhnear equation 
which is the form 

(8) 



"' = 7- 



where g {z, t) are complex functions and / {z, t) is a real function. Substituting ([8]) into 
we obtain 



-/j (^Dt-^Dl]g.f-4Sg 



Dlf -f + ^gg 



S"" \ASJ 

where Dt and Dj. are the Hirota bilinear operators which are defined as 

Then Eq. ([9]) can be decoupled as two equations 

Gig-f = 0, 

3 

G2f ■ f = -^99, 

where the Hirota bilinear operators Gi and G2 are given by 

2JS 



2 f 



0, 



(9) 



(10) 



z=z',t=t' 



(11) 



h 



-Dt - A, 



G2 = Dl+^'' 



(12) 



where A is constant to be determined. Now the Eq. ffTTl) has made the Eq. ([7]) to the 
normal procedure of Hirota method for getting the exact soliton solutions. In the following 
by making a series of suitable assumption for the expression of g and /, the exact one- and 
two-soliton solution can be obtained analytically. To this purpose we suppose firstly that 



^ = ^0 (1 + X9i) , / = 1 + X/: 



1) 



(13) 



in order to construct a dark one-soliton solution for the system ([7]), where x is an arbitrary 
parameter which is absorbed in expressing the soliton solutions in the following sections. 
Substituting f|T3l) into ffTTl) . and collecting the coefficients of x", we obtain 



Gigo ■ 1 



0, 



ASJ ASJ 



S 



9o9o, 



(14) 
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where the first equation in Eq. f ll4p can be expressed by the definition fllUp from which in 
is easy to find the solution satisfying (H^ as 

5(0 = Tiexp (i$i) , (15) 

where $1 has of form 

$1 = ][^hz -{\-ll)t- (16) 

in which li, $5 ^-^'s real constants and ri are complex constant. With the restriction, i.e., 
the second equation in Eq. (fT^ the expression of ri is 

\T^f = ^{lB, + 8SJ-hX) (17) 

which imply the condition '-/Bz + 8SJ — hX > 0, i.e., the minimum external field value for 
the existence of soliton solution in Eq. (jZj), i?^ = (^A — 8SJ) /•y. The coefficient of x in Eq. 
(fTTD leads to 



£i(l-/i + (7i-l)=0, 

's 



G2(/i-l + l-/i) = -||n|'(^?i+^i), (18) 



where £1 is defined as 



£1 = iD, - - z2h^^Dz. (19) 



By expanding Eq. (|T8|) with the definition (|T0|) in detail, one can easily check that equations 
f lTSjl admit the following solutions 

h = exp^i, (20) 

where 

in which Pi, fiiand are real constants and Zg is a complex constant which is given by 



Z 



(21) 
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From this equation it can easily be seen that | | ^ = 1 • Then using Eq. fll3p , after absorbing 



X, the dark one-sohton solution of spinor BECs in an optical lattice can be derived as 



•"■ = 7 = T 



:i + z,) - (1 



Zg) tanh ^ 



exp [i^i] 



(22) 



The solution in Eq. ( 122|) describes a magnetic dark soliton solution characterized by four 
real parameters: velocity of envelope motion fii/Pi, phase $i, coordinate of the center of 
the solitary wave ^i/Pi and initial phase ^^^^$5/^1- This solution is similar to that of the 
Heisenberg spin chain with an external field where the dipolar coupling is typically several 
orders of magnitude weaker than the exchange coupling and would correspond to Curie 
temperatures much below the observed ones. Hence its contribution to the spin wave and 
magnetic soliton can be neglected in practice. However, for the spinor BECs in the optical 
lattice the exchange interaction is absent and the individual spin magnets are coupled by the 
magnetic dipole-dipole interactions. Due to the large number of atoms at each lattice site, 
these site to site interactions, despite the large distance between sites, explain the natural 
existence of magnetic soliton which agrees with the results in Refs. . With the help of 
Eq. (ITSi) we have 

(23) 



By combine the above presentation with the Eq. (|T7|) we obtained the following relation for 



5. 



1 
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2/iv^+^=-^) -8SJ + hX 
V2JSP1J 



which shows that the external field Bz and the dipole-dipole interaction can affect the 
velocities, size and shape frequency of the magnetic soliton. It offer an useful scheme to 
control soliton in practical experiment by adjusting of the magnetic field. 



IV. DARK SOLITON INTERACTION 

In this section we will give the analytical expression of two soliton solutions of Eq. ([7]). 
The properties is discussed in detail as well. To this purpose for constructing the dark 
two-soliton solutions of Eq. ([7]) we now assume that 

g = go{l + X9i + X^92) , / = 1 + x/i + x'f2, (24) 
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where go and Hq are obtained here as in Eq. (IT^ . By employing the same procedure before 
we obtain the following set of equations from Eq. f lTTj) . corresponding to the different powers 
of X as follows 

(i) for the coefficient of x 

G2(/i-l + l-/i) = -||ri|'((7i+^i), (25) 



where £i is same with Eq. (HM . 
(ii) for the coefficient of x^ 



Gi {go ■ /2 + goQi ■ fi + 9092 • 1) = 0, 
£i{l-f2 + 9i-fi + 92-l) = 0, 



G2 (1 ■ /2 + /2 ■ 1 + /i ■ /i) = -I [kil' toi + 92 + 92)] ■ (26) 



S 

(iii) for the coefficient of x^ 



Gi {9o9i ■ f2 + 9o92 ■ fi) = 0, 

£1 (91 ■ f2+ 92- fl) = 0, 

G2 {fl ■ /2 + /2 ■ fl) = -| \ri\' {91-92 + 92-gi) . (27) 
(iii) for the coefficient of x^ 

Gi {9092 ■ f2) = 0, 

£1 {92 ■ f2) = 0, 

G2{f2-f2) = -^\rif92g2- (28) 

Repeating the same normal procedure in obtaining the one-soliton solution one can easily 
get the solutions of the above set of equations as 

gi = Ziexp^i + Z2exp^2> 
g2 = A12Z1Z2 exp (^1 + C2) , 
fl = exp ^1 + exp ^2, 

/2 = Ai2exp(ei+^2), (29) 
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where 



z 



3 



Pj + zJ6\T,f/S-pf 



in which j = 1,2. Noting that 



_ \ ■ 2JS 

s h ' 

by which we can rewrite Zj and fi-,-, j = 1, 2, as 



1/2 



-'3 



'-fP,+ZKi^ 



= \n,-2hJ^\P,, 



2JS 



So from Eq. (125|) to Eq. (I28i) the expression of A12 is obtained in the form 

{P,-P,fm^f + {n,-K,f 



A 
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{P, + P.f{^Y + {K,-K,f 

Now using Eq. ([8]), ( 124|) and ( 129|) . the dark two-sohtons can be found exphcitly as follows 

/ 3 l + Ziexp^i + Z2exp^2 + ^i2^i^2exp(^i+^2) /gn^ 

^t«;o = 7 = n exp ?$i — — (c ^ c \ ' 30 

/ 1 + exp ^1 + exp ^2 + ^12 exp (^^ + ^2) 



where all the parameters have been obtained in the previous section. The solution (130|1 
describes a general elastic scattering process of two solitary waves with different center 
velocities Vti/Pi and ^2/^2 and the same phase $1 . Before collision, they move towards 
each other, one with velocity Vti/Pi and the other with Q2/P2- In order to understand the 
nature of two-soliton interaction, we analyze the asymptotic behavior of two-soliton solution 
(130|) . Asymptotically, the two-soliton waves ( 130|) can be written as a combination of two 
one-soliton waves fl22l) . The asymptotic form of two-soliton solution in limits t — )■ —00 and 
t — )■ 00 is similar to that of the one-soliton solution (|22|) . 

(i) Before collision (limit t — )■ —00). 

(a) Soliton 1 (^^ ^ 0, ^2 ~^ —00). 



i^two n exp(i$i) 



:i + Zi)-(l-Zi)tanh 2 



(31) 
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(b) Soliton 2 ( ^2 ~ 0, ^1 ^ oo). 



;i + Z2) - (1 - Z2) tanh 



^2 + ^0 



(32) 



where the initial coordinate of the center of the sohtary wave is removed by 60 = lnAi2. 
(ii) After colhsion (hmit t — ?■ 00). 



(a) Sohton 1 {^^ ^ 0, ^2 

i'two ^2^1 exp {i^i) 

(b) Sohton 2 ( ^2 ~ 0, ei ^ -00). 



'1 + Zi) - (1 - Zi)tanh 



^1 + ^0 



(33) 



1 - Z2) tanh y 



(34) 



Analysis reveals that there is no amplitude exchange for soliton 1 and soliton 2 during 



colhsion because of \Zj\ = 1, j = 1,2. However, from Eqs. fl3T|) to f lMjl one can see 
that there is a phase exchange Sq/2 for soliton 1 and soliton 2 during collision. It shows 
that the information held in each soliton will almost not be disturbed by each other in 
soliton propagation. These properties may have potential application in future quantum 
communication. 



V. CONCLUSION 

Dark soliton dynamics of spinor BECs in an optical lattice is studied in terms of a 
nonlinear schrodinger equation by means of Hirota method. Exact soliton solutions are 
obtained analytically and the elastic collision of two solitons is demonstrated. It should be 
interesting to discuss how to create and to detect such magnetic soliton in experiment. Using 
Landau-Zener rf-sweeps at high fields (30 G) ^ a condensate was prepared in the hyperfine 
state 1/ = lyTTLf = 0) of sodium, i.e. the ground state of the spinor BECs. Then the atoms 
of the ground state can be excited to the hyperfine state |/ = l,mf = ±1) by laser light 
experimentally. Therefore the excited state of the spinor BECs, i.e. the magnetic soliton 



can be created. As the same discussion before [12| the spatial-temporal spin variations in 
the soliton state are significant. This makes it possible to make a direct detection of the 
magnetic soliton of spinor BECs. By counting the difference numbers of the population 
between the spin +1 and —1 Zeeman sublevel, the average of spin component < 5*^ > is 
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measured directly. Transverse components can be measured by use of a short magnetic 
pulse to rotate the transverse spin component to the longitudinal direction. Any optical 
or magnetic method which can excite the internal transitions between the atomic Zeeman 
sublevels can be used for this purpose. In current experiments in optical lattices, the lattice 
number is in the range of 10-100, and each lattice site can accommodate a few thousand 
atoms. This leads to a requirement for the frequency measurement precision of about 10-100 
kHz. This is achievable with current techniques. 

The magnetic sohton of spinor BECs in an optical lattice is mainly caused by the magnetic 
and the light-induced dipole-dipole interactions between different lattice sites. Since these 
long-range interactions arc highly controllable the spinor BECs in optical lattice which is an 
exceedingly clean system can serve as a test ground to study the static and dynamic aspects 
of soliton excitations. 
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